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We study the possibility to extract model independent information about the dynamics of the 
universe by using Cosmography. We intend to explore it systematically, to learn about its limi- 
tations and its real possibilities. Here we are sticking to the series expansion approach on which 
Cosmography is based. We apply it to different data sets: Supernovae Type la (SNela), Hubble 
parameter extracted from differential galaxy ages, Gamma Ray Bursts (GRBs) and the Baryon 
Acoustic Oscillations (BAO) data. We go beyond past results in the literature extending the series 
expansion up to the fourth order in the scale factor, which implies the analysis of the deceleration, 
go, the jerk, jo and the snap, so. We use the Markov Chain Monte Carlo Method (MCMC) to 
analyze the data statistically. We also try to relate direct results from Cosmography to dark energy 
(DE) dynamical models parameterized by the Chevalier-Polarski-Linder (CPL) model, extracting 
clues about the matter content and the dark energy parameters. The main results are: a. even if 
relying on a mathematical approximate assumption such as the scale factor series expansion in terms 
of time, cosmography can be extremely useful in assessing dynamical properties of the Universe; b. 
the deceleration parameter clearly confirms the present acceleration phase; c. the MCMC method 
can help giving narrower constraints in parameter estimation, in particular for higher order cosmo- 
graphic parameters (the jerk and the snap), with respect to the literature; d. both the estimation of 
the jerk and the DE parameters, reflect the possibility of a deviation from the ACDM cosmological 
model. 

PACS numbers: 04.50.-h, 98.80.-k, 98.80.Es 



I. INTRODUCTION 

Cosmography is a branch of cosmology which has not 
been explored so much yet in spite of its promising 
prospects towards a deeper understanding of the accel- 
eration of the universe. Thus, there remain caveats to 
be filled and misunderstandings to be rectified before one 
can take full advantage of the possibilities offered by Cos- 
mography. 

This approach proceeds by making minimal dynami- 
cal assumptions, namely, it docs not assume any partic- 
ular form of the Friedmann equations and only relies on 
the assumption that the spacetime geometry is well de- 
scribed on large scales by the Robertson- Walker metric. 
This makes cosmography fully model and setting inde- 
pendent, in the sense that it is only at the late stages of 
interpretation of the results offered by this approach that 
one has to care for the theoretical framework in which the 
conclusions fit best (dark energy, dark gravity, modified 
gravity, etc.). 

This is clearly important, for instance, when one 
attempts at comparing the ACDM model (or any dark 
energy model) with alternative theories of gravity (such 
as f{R) gravity). In the traditional approach a parame- 
terized model is assumed at the beginning, and then, by 
contrasting it against the data, it is possible to check its 
viability and to constrain its characterizing parameters. 
But this approach is clearly a model dependent one 
so that some doubts always remain on the validity of 



the constraints on derived quantities as the present day 
values of the deceleration parameter and the age of the 
universe, just to mention a couple of them. Conversely, 
the constraints that one can infer from cosmography are 
more universal since they are completely model-free. If 
well understood and analyzed in full details, weighting its 
pros and cons, results from Cosmography can be rightly 
considered as milestones in the study of properties of 
universe dynamics, which any theoretical model has to 
consider and satisfy, and can help pointing in the right di- 
rection in connection with the properties of the Universe. 

Of course, everything is not so clearly cut. The fun- 
damental rule of cosmography is based on the series ex- 
pansion of the scale factor, a(t). From it one obtains the 
cosmographic parameters usually referred to as the Hub- 
ble (H), deceleration (q), jerk (J), snap (s) and lerk (I) 
parameters, which we define below. However, the series 
of a(t) places difficulties almost from the very beginning. 
As discussed in [![, there are two main problems when 
working with series expansions: the convergence and the 
truncation of the series. Those authors show it is possible 
to solve the convergence problem defining a new rcdshift 
variable, the so called y-redshift: 



For a series expansion in the classical z-redshift the con- 
vergence radius is equal to 1 . This represents a drawback 
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when one wants to extend the application of cosmogra- 
phy to redshift z > 1, which is not such a big past ex- 
tension, if we consider that current supernovae data go 
up to z = 1.4, or that the CMB analysis involves the 
decoupling redshift z ps 1089. It would be pleasant to be 
able to work correctly at least up to the highest redshift 
demanded by the observational data set one wants to use. 

The y-redshift could potentially solve this problem be- 
cause the z-interval [0, oo] corresponds to the y- interval 
[0, 1], so that we are mainly inside the convergence in- 
terval of the series, even for CMB data (z = 1089 — > 
y = 0.999). So, in principle, we could extend the se- 
ries up to the decoupling redshift value, and one could 
place CMB related constraints within the cosmographic 
approach. Unfortunately, some problems arise and we 
will discuss them below. Leaving aside this limitation, 
the possibilities offered by the y-redshift can be consid- 
ered as a step in the right direction, because, obviously, 
defining a new redshift does not change the physical con- 
tent of the data. 

But the problem of the truncation of the series is not 
solved; of course, we expect that for low redshifts a low 
order series expansion will work well, while for higher 
redshifts it is likely that higher order expansions of any 
physical quantity of relevance will be required. But, is 
this true? And what is the good order for any redshift? 
Is there an even approximate relation which correlates 
redshift with series order? 

This paper represents an improvement in this direc- 
tion. As considering a higher order of expansion enforces 
taking into account more free parameters, it seems a pri- 
ori it will be difficult to obtain tight constraints on them 
if the discriminating power of the data sets used is weak 
or not fully exploited statistically. Despite this, we suc- 
cessfully put bounds on the parameters of the third order 
series in a statistically consistent way, and so make a con- 
siderable improvement with respect to @, where for the 
same redshift extent the authors only managed to put 
constraints on the coefficients of a second order series. 
This is one of the major contributions of our paper. 

In addition, there is another matter on which this pa- 
per enhances the worth of the cosmographic approach: 
the values and the errors of the cosmographic parameters. 
Estimations of that sort can be found in the literature, 
not only in Q, but also in 1 and 0. While the former 
used supernovae data, in related works like [B-Q one can 
find theoretical estimations (in the context of f(R) theo- 
ries in both metric and Palatini approaches) which could 
be related to estimations from observational data. Some 
recent contributions to this topic are the works [T3, EH , 
where gamma ray burst data have been used. Here we 
analyze at the same time data from SNela and GRBs, 
and combine them with a sample of Hubble parameters 
values extracted by means of the differential galaxy ages, 
and BAO data. The use of the Hubble parameter data, 
in particular, introduces additional advantages we will 
refer to later. 

The main problem with earlier contributions to the 



topic is that the error bars obtained were considerably 
large. Due to this fact, the statistical conclusions to 
be inferred are rather weak and, in consequence, they 
hardly provide valid priors for future analysis. For exam- 
ple, from [| we have q = -0.90 ± 0.65, j = 2.7 ± 6.7, 
.s = 36.5 ± 52.9, l = 142.7 ± 320; some of the rela- 
tive errors are as high as about a 200% (as usual, the 
subindcx denotes the z = values of the parameters). 
Now, we can expect that errors on the coefficients will 
get bigger when higher order expansions are considered. 
As those parameters are correlated among them, errors 
in the low order series coefficients propagate to the addi- 
tional coefficients included in the higher order series. If 
these errors turn out to be big, the positive aspects of 
cosmography will vanish. For instance, if the error on qo 
allows for a positive value of the deceleration parameter 
in the 3cr confidence interval, it will not be so obvious 
that the Universe is accelerating (see again Q). More- 
over, the ACDM model enforces jo = 1 (this can be easily 
derived from the definition of jo , which is to be found in 
§ piip ): but if the error bar is as big as reported above 
(from [J]), it will not be possible to confront this model 
with competing ones. 

Given those problems we have just discussed, an im- 
portant question arises: is it not possible to give narrower 
statistical constraints? To this end we apply a Markov 
Chain Monte Carlo (MCMC) method: it allows us ob- 
taining marginalized likelihoods on the series coefficients 
from which we infer rather tight constraints on those pa- 
rameters. The reason for this considerable improvement 
with respect to earlier works is that in our code we have 
implemented several controls which give us power over 
any physical requirement we expect from our theoretical 
apparatus. Actually, setting restrictions on the Hubble 
parameter is possible because we use data related to this 
quantity that reveal to allow a considerable improvement 
in the quality of constraints. 

We have organized the paper as follows: in § ((TT| and 
Appendix) we define the cosmographic parameters and 
give all the relations needed; in § we describe some 
caveats that have to be taken in mind when working with 
Cosmography; in § (|IV[) we describe the observational 
data used for the analysis; in § (JVj and § (|VI[) we present 
our main results and discuss their meaning and conse- 
quences; finally, in § (|VII[) we summarize. 



II. COSMOGRAPHIC PARAMETERS 

The key rule in cosmography is the Taylor series expan- 
sion of the scale factor with respect to the cosmic time. 
To this aim, it is convenient to introduce the following 
functions: 



1 d 2 a 1 
q{t) = —al*lP' 
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1 d 3 a 1 
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They are usually referred to as the Hubble, deceleration, 
jerk, snap and lerk parameters, respectively (see [l2| for 
an historical account of these names). Using these defi- 
nitions it is easy to write the fifth order Taylor expansion 
of the scale factor: 



a (t ) 



= 1 + H (t t ) |i? 2 (t - t ) 2 + |i? 3 (t - t ) 3 
+ ^H 4 (t - t ) 4 + l ^H 5 (t to) 5 + 0[(t to) 6 ] (3) 



with to being the current age of the universe. Note that 
Eq. is also the fifth order expansion of (1 + z) -1 , as 
the definition for the redshift z is z := a(t )/a(t) — 1. 
This is what we need for developing the cosmographic 
apparatus; for sake of clearness we remind all the detailed 
calculations to the appendix section, while giving here 
only the main results. 



A. Luminosity distance series 

Since we are going to use SNela and GRBs data, 
it will be useful to give the Taylor series of the ex- 
pansion of the distance modulus, which is the quantity 
about which those observational data typically inform. 
The final expression for the distance modulus based on 
the Hubble free luminosity distance (Eqs. (|A5|) - (|A6|) ^ . 
fi(z) = 51og 10 di,{z) + fiQ, is the following: 



B. Hubble parameter series 

The definition of the luminosity distance given by 
Eq. (|A3|) . can be presented in this other way: 



D L {z) = c (l + z) 



(5) 



It is interesting to study the possibility of obtaining the 
same final expression for the distance defined starting 
from a Taylor series expansion of the Hubble parameter 
instead of the scale factor, namely: 



H(z) = H Q + ^- 
dz 



1 d 2 H 
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This series expansion will be also an important rule in 
the definition of our Markov chains algorithm given the 
observational data we use. To compute all the terms of 
this series we have to keep in mind the derivation rule 
(for a clearer notation we will suppress the redshift de- 
pendence of the Hubble parameter setting H = H(z)): 



d (1 + )H ^ 
dt dz 



(7) 



It is an easy but cumbersome task to obtain from the 
latter the higher-order time-derivatives of the Hubble pa- 
rameter (they are given in § ([U]). Then, from the defini- 
tions of the cosmographic parameters, Eq. (j2J, it is easy 
to demonstrate the following relations: 



H = -H 2 (l + q), 



H = H 3 (j + 3q + 2). 



(8) 
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fx(z) 



with 



log 10 



(log z + M 4 z + M 2 z 2 + M 3 z 3 M 4 z 4 ) 
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M 1 = --[-l + q ] , 

M2 = -^ [7-10«)-9gg+4jo] , 

M 3 = i [5 - 9q Q - I6q 2 - 10q 3 + 7j + 8q j + s ] 



M 4 



1 



[-469 + 1004 go + 2654c? 2 . + 3300^ + 1575g, 



2880 
200jj 

324s - 24Z 



200j 2 - 1148jo - 2620 9o jo - 1800ggj - 300^0 



Of course we have also derived the same relations for the 
y-redshift, but we relegate them to the appendix. 



d 3 H 



dt 3 
d 4 H 



= H 4 [s - 4j - 3q(q + 4) - 6] 



dt 4 



= H 5 [I - 5s + 10(q + 2)j 
+30(g + 2) g + 24] 



(10) 



(11) 



If we convert time-derivatives into derivatives with re- 
spect to redshift using Eqs. - (jCTJ) - |C2j - JC3} wc 
have: 



dH H M 

-f- = T— (1 + 9) 

dz l + z 



d 3 H 



dz 3 (1 + 2) 3 



d 2 H H 
llz 2 ~ (l + z) 2 

H 



(12) 
(13) 



(3q 2 + 3<? 3 - 4qj - 3j - s) , (14) 
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d 4 H H 
~dz T ~ (1 + z 

+ 7gs + 12j - Aj 2 + 8s + I) 



- {-I2q 2 - 2Aq 3 - 15g 4 + 32qj + 25q 2 j 

(15) 



Exactly the same luminosity distance formula as given 
in Eq. (|A5|) can be recovered after a series of steps: after 
having evaluated Eqs. (JT2|) - (|15|) at z = (namely to), 
we must insert them in Eq. (j6|), and then insert this one 
in Eq. ([SJ. Then we must Taylor expand the integration 
function, integrate it, and we will have reached the final 
sought result and the coincidence will be complete. 

This is all, so far, in what the basic setup is concerned, 
in the next section we move on to applications of these 
definitions. 



III. SERIES ORDERS COMPARISON 

As we have said in the introductory section, the first 
question we would try to answer is whether there is a re- 
lation between the highest expansion order in the Taylor 
series and the redshift range where this series can be ap- 
plied. What one could expect is the rough rule that we 
would need a series expansion truncated at higher orders 
when increasing up the redshift range. This would be a 
mathematical requirement, of course. But what about 
physics, our main interest, and the fit to the Hubble 
SNela diagram; namely, the reproduction of the luminos- 
ity distance or its corresponding distance modulus? We 
have found that the answer seems to be not so obvious. 

Let us start by writing some useful relations for suc- 
cessive results. To make the discussion more specific we 
will use the Chevallier - Polarski - Linder (CPL) ([3 , (3 ) 
paramctrization: 



w = wo + w a (l - a) = w + w a z(l + z)~ 



(16) 



so that, in a spatially flat universe filled with cold dark 
matter and dark energy, the dimensionless Hubble pa- 
rameter E(z) = H/Hq reads: 



e 2 (z) = n m (i+ z ) 3 +n x (i+z) 



(17) 



with fix = 1 — because of the flatness assumption. In 
order to determine the cosmographic parameters for such 
a model, we avoid integrating H(z) to get a(t) by noting 
that d/dt = — (1 + z)H(z)d/dz. We can use such a re- 
lation to evaluate (H,H,d 3 H/dt 3 ,d 4 H/dt 4 ) and convert 
Eqs. flSJ - (fTTj) to redshift derivatives, i.e.: 



3{z) 



(l + z) dH 2 (z) 



2H 2 {z) d. 
(l + z) 2 d 2 H 2 {z) 
2H 2 (z) dz 2 



-1, 



(l + z) dH 2 (z) 
H 2 (z) dz 



1.(18) 



For the snap and the lerk the expressions are much longer 
and involve respectively the third and the fourth order 
derivatives of H 2 (z). We underline that these expres- 
sions for cosmographic parameters are exact. Given an 



analytical expression for H(z), we can evaluate them at 
z = 0, and solve them with respect to the parameters of 
interest. 

Some algebra finally gives: 



1 3. . 



(19) 



jo = l + -(l-n m )[3wo(l + w )+w a ] , (20) 



7 33 n o ^ 
so = -- - —(1 - il m )w a 

- j(i - Jim) [9 + (7 - n m )w a ] wo - |(i - n m ) 

27 

x (i6 - m m )w 2 —{i - n m )(3 - n m )w 3 Q , (21) 



h = y + [213 + (7 - a m )w a ] w a 



[489 + 9(82 - 2ltt m )w a } w 



+ gCl-fim) 

27 



67- 21fL 



-(23-110 m K 



— (i-n m )(47-24n w K 

81 

— (l-Q m )(3-2fi, N 



mJuJQ 



(22) 



From Eq. (fTT)) we can derive the exact analytical ex- 
pression for the Hubble free luminosity distance and the 
distance modulus. Then, we will compare it with its ex- 
pression as a cosmographic series, i.e. a series depending 
on the cosmographic parameters. We have considered 
the cosmographic series with two (qo,jo), three (qo , jo , s o ) 
and four parameters (go, jo, so, lo)- Of course, from any 
CPL scenario we will be able to derive the corresponding 
set of cosmographic parameters using Eqs. (fT9"f - ([22]) . 

We have tested two different toy models: a ACDM 
model (£l m = 0.3 , wq = —l,w a = 0); and a dynamical 
dark energy model with Cl m = 0.245 , wo — —0.93 , w a = 
—0.41 coming from (l5| . In Figs. ((T|) - ^ we show the 
comparison in the redshift range < z < 2, which is 
the range where cosmography has been mainly applied in 
past works, regarding to the maximum available redshift 
of the SNela surveys used in the literature. 

As we can see in Fig. (JlJ, it is not so obvious that 
the series expanded at higher orders can describe well 
the underlying Hubble parameter or distance modulus 
when applying cosmography out of the convergence ra- 
dius. There we show the relative residuals with respect to 
the exact CPL functions (i.e. the Hubble parameter and 
the distance modulus) coming from the different order 
series in the ACDM case. As we can see, inside the red- 
shift convergence radius (z < 1) the four-parameter series 
effectively gives a better approximation to the Hubble 
parameter and the distance modulus than other orders 
expansions. 
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Figure 1. (Top.) Residuals between the CPL exact Hub- 
ble function and different-order series expansions of the same 
quantity in the ACDM case. The redshift range is < z < 2. 
(Bottom.) The same as above picture but for distance mod- 
ulus. One might except that the value of the denominator, 
(J.ACDM, depends on the nuisance parameter fj,o, i.e. on the 
value of the Hubble constant, Ho. But we have checked that 
changing Ho in the range [0.65,0.75], which is large enough 
to contain all the acceptable physical values of Ho, does not 
produce any sensible change in the figures. 



But it is also clear that the rough rule (higher redshift, 
higher order series) is somewhat not followed when out 
of the convergence interval, namely, above 2=1. In this 
case residuals from the four-parameter expansion grow 
with redshift faster than the ones coming from the three 
and two-parameters ones. 

The same analysis made with the dynamical model 
lead to similar final conclusions and adds further inter- 
esting caveats which have to be kept in mind when work- 
ing with cosmography and using its results. Looking at 
Fig. ^ we can see some consequences from changing the 
cosmological model: 

• the inversion in goodness between the order series 
expansions even inside the convergence radius. In 
the H(z) case, we can see that in the ACDM model, 
for z > 1, the two and the tree-parameters expan- 
sions were quite similar, with a slight preference 
for the latter one. In the dynamical dark energy 
model things go different, and the two-parameters 
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Figure 2. (Top.) Residuals between the CPL exact Hub- 
ble function and different-order series expansions of the same 
quantity in the model from [l5l |. The redshift range is 
< z < 2. (Bottom.) The same as above picture but for 
distance modulus. 



expansion is strongly preferred; 

• different magnitudes for residuals, notably in the 
Hubble function case, slightly in the distance mod- 
ulus. 

The first point leads us to think about the real ef- 
fectiveness of cosmography approach as a "model inde- 
pendent" one. Cosmography is model independent, of 
course, starting from its basic assumptions, and the in- 
formation pieces one can extract from it are model in- 
dependent too; but data are tracers of a cosmological 
model. 

From the previous considerations it seems that there is 
an implicit relation between the series expansion and the 
underlying unknown cosmological model which we want 
to read out from data. If cosmography, by its order series 
expansion, is sensible enough to detect and discriminate 
a model against another one, then the choice of working 
with two, three or four parameters can be critical. 

The second consideration is important for defining the 
main lesson of this section. All the analysis we have 
done until now could be useless without relating it to the 
experimental possibilities we have nowadays. Namely: 
are our present data and their related errors able to make 
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us distinguish between a two, three, or four-parameter 
expansion for the measured physical quantities we are 
working with? 

To answer this question we need a look at the data we 
are going to use. If we consider the data from [l6|, we 
see that the relative errors on H(z) range in the inter- 
val [0.10,0.62], which is clearly larger than the residual 
we have plotted in Figs. dTJ - ((2j) . So any possibility to 
discriminate what the right expansion order is in cosmog- 
raphy is completely out of question. Things are slightly 
better with the SNela data set, which shows relative er- 
rors in a range [0.002, 0.025]. But even now results have 
to be taken with a pinch of salt. For example, in the 
rcdshift range z < 1, SNela relative errors are < 0.010, 
which means that they are of the same order of the resid- 
uals we have found. In this case, we can say we are in a 
"border-line" situation. 

Last but not least, we have done the same analysis 
working with the y-redshift. In this case we have found 
out that the higher order series give a better approxi- 
mation to the exact relation than the lower order ones, 
which is quite expected, because we are always well in- 
side the convergence radius, i.e. < y < 1. On the 
contrary, the most important thing that emerges, is that 
the y-redshift can be useful only up to a certain redshift. 
Above a certain limit it cannot be used to constrain the 
Hubble parameter or distance modulus, for some intrin- 
sic mathematical properties coming from its definition. 
As it is showed in the Fig. ([3]), above y ps 0.4 0.5, 
there is a clear deviation between the series expansions 
and the exact expression of our physical quantities, larger 
than observational errors. And the y-redshift series are 
unable to follow the right trend. 

For all these reasons, we have chosen to work in the 
reminder with second and third order series expansions, 
giving at most constrains to the deceleration, the jerk and 
the snap parameters. We have seen that it is impossible 
to extend the analysis to four parameters, the lerk being 
completely undetermined. 



IV. OBSERVATIONAL TESTS 

A. Hubble parameter 

Recently in [l6| an update of the Hubble parameter 
H(z) data extracted from differential ages of passively 
evolving galaxies previously published in [l?} was pre- 
sented. In |l8| this method was analyzed and some 
caveats about its use are exposed; anyway, data in (l6j 
are considered greatly improved with respect of previous 
ones derived in [l7j . Constraining the background evo- 
lution of the universe using these data is interesting for 
several reasons. Firstly, they can be used together with 
other cosmological tests in order to get useful consistency 
checks or tighter constraints on models. Secondly, it is 
well known that detailed and reliable information in the 
behavior of H(z) is lost when using the luminosity dis- 




Figure 3. (Top.) Hubble parameter in the ACDM model. 
The y-redshift range is < y < 1. The solid black line is 
the exact analytical expression; the dashed line is the fourth 
order series; the dot-dashed line is the third order series; and 
the dotted line is the second order series. The vertical dotted 
line is the maximum redshift from SNela sample; the vertical 
dotdashed line is the maximum redshift from Hubble data; 
the vertical dashed line is the maximum redshift from GRBs 
sample. (Bottom.) The same as above picture but for distance 
modulus. 



tance or the correlated modulus distance. As a result of 
the integration of this function, the fine details cannot be 
analyzed and no information can be derived (if there is 
any information, probing or refuting the previous ones). 
For this reason, we have chosen to work with this data 
set: any refinement in the Hubble parameter could give 
important details. 

The Hubble parameter dependence on the differential 
age of the Universe in terms of rcdshift is given by 



H(z) = 



1 dz 
1 + zdt' 



(23) 



Thus, H(z) can be determined from measurements of 
dt/dz. As reported in [l6j|, [l7j], [H and Q, values of 
dt/dz can be computed using absolute ages of passively 
evolving galaxies. 

The galaxy spectral data used by [T(| come from 
observations of bright cluster galaxies done with the 
Keck/LRIS instrument (see 2l[ for a detailed description 
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of the observations, reductions and the catalog of all the 
measured redshifts). The purposely planned Keck-survey 
observations have been extended with other datasets: 
SDSS improvements in calibration available in the Public 
Data Release 6 (DR6) have been applied to data in [20| : 
the SPICES infrared-selected galaxies sample in [22f: and 
the VVDIS survey by the VLT/ESO telescope ing3|. 

The authors of these references bin together galaxies 
with a redshift separation which is small enough so that 
the galaxies in the bin have roughly the same age; then, 
they calculate age differences between bins which have 
a small age difference which is at the same time larger 
than the error in the age itself (|16|). The outcome of 
this process is a set of 12 values of the Hubble parameter 
versus redshift. A particularly nice feature of this test 
is that the sensitivity of differential ages to systematic 
error is lower than in the case of absolute ages (pijy 

Observed values of H(z) can be used to estimate cos- 
mographic parameters by minimizing the quantity 



^ (H(z ]] {6 l )}-H obs (z ] ))' 



°h(*j) 



(24) 



where H = 100 h will be fixed as h = 0.742 (|25j), while 
6i is the vector of model parameters, namely in our case 

0i = (go, jo, so)- 

To minimize the x 2 function we will use Markov Chain 
Monte Carlo (MCMC) methods (fully described in [H, 
[27j , [28| and references therein) testing their convergence 
with the method developed and fully described in |29| . 



B. Supernovae 

We use the most updated SNela sample we have now, 
the Union2 sample described in [3(| • The Union2 SNela 
compilation is the result of a new low-redshift nearby- 
Hubblc-flow SNela and new analysis procedures to work 
with several heterogeneous SNela compilations. It in- 
cludes the Union data set from [3l| with six SNela first 
presented in [3(|, with SNela from f32|. the low-z and 
the intermediate-z data from [33| and |34| respectively. 
After various selection cuts were applied in order to cre- 
ate a homogeneous and high signal-to- noise data set, we 
have final 557 SNela events distributed over the redshift 
interval 0.015 < z < 1.4. 

The statistical analysis of Union2 SNela sample rest 
on the definition of the distance modulus, 



fi(Zj) = 5log w (d L (Zj, {9i})) + fio 



(25) 



where d^iz^ {0i}) is the Hubble free luminosity distance, 
Eq. (|A5[) , expressed as a series and depending on the 
cosmographic parameters, 9i = (goi joi s o)- The best fits 
were obtained by minimizing the quantity 



XsnOo, 



E 7 (M^MO, ~ Vobs{Zj)Y 

a 2 - 

3 = 1 M 



(26) 



where the a 2 „■ are the measurement variances. The nui- 
sance parameter /j,o encodes the Hubble constant and 
the absolute magnitude M, and has to be marginalized 
over. Giving the heterogeneous origin of Union data set, 
and the procedures described in 3l| for reducing data, 
we have worked with an alternative version of Eq. (|26[) , 
which consists in minimizing the quantity 



_1 

C3 



with respect to the other parameters. Here 



Cl 



557 

£ 



OO^SMo = 0, {6»j)} - Hobs(zj))' 2 



(27) 



(28) 



557 



c 2 - -ji > y^i 

= l 



C3 



557 
3=1 



(30) 



It is trivial to see that Xsn ^ s J us ^ a version of Xsn-> 
minimized with respect to /i - To that end it suffices to 
notice that 



Xsn(mo, = ci - 2c 2 /i + c 3 /Xq 



(31) 



which clearly becomes minimum for fio = C2/C3, and so 
we can see Xsn = X§n(mo = 0,{^i})- Furthermore, one 
can check that the difference between Xsn an d Xsn ^ s 
negligible. 



C. Gamma Ray Bursts 

Working and interpreting results from a GRBs analysis 
is not an easy task, being the errors on their observable 
quantities much larger than those ones for SNela, and 
being their source mechanism still not well understood. 
For this reason choosing a good GRBs sample is crucial; 
we have chosen to work with the sample described in 
[35l] . There the authors perform a new calibration pro- 
cedure on the widely used GRBs sample from [36| which 
perfectly matches our requirements for using them with 
cosmography. The possibility of a joint analysis of SNela 
and GRBs is strictly related to the building of a Hubble 
diagram for GRBs too, being this extremely difficult be- 
cause of GRBs are not standard candles as SNela. To 
create an Hubble diagram for GRBss, one has to look 
for a correlation between a distance dependent quantity 
and a directly observable property. Starting from some 
of the many correlations that have been suggested in the 
last years, [36| created a sample of 69 GRBs in the red- 
shift range 0.17 < z < 6.6 whose Hubble diagram is well 
settled. 
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In [351 ] the authors have updated such a sample in many 
aspects, the main one being the test of a new method 
for the calibration of GRBs based one the assumption of 
none a priori cosmological model. Such a model inde- 
pendent calibration is built on the idea that SNela and 
GRBs at the same redshift should exhibit the same dis- 
tance modulus. In this way, interpolating the SNela Hub- 
ble diagram gives the value of n(z) for a sub-sample of 
GRBs which lies in the same redshift range. This sub- 
sample can be finally used for calibrating the well know 
GRBs correlations and, assuming that this calibration is 
redshift independent, it can be extended to high redsfhit 
GRBs. With this procedure, in [35| the authors were 
able to convert the (36j sample to a new one, with the 
same number of objects but with SNela-calibrated GRBs 
distance modulus. These data can then be used for min- 
imizing the corresponding x 2 : 



3=1 



a 



(32) 



D. Baryonic Acoustic Oscillations 

In [3?! the authors analyze the clustering of galax- 
ies within the spectroscopic Sloan Digital Sky Survey 
(SDSS) Data Release 7 (DR7) galaxy sample, including 
both the Luminous Red Galaxy (LRG) and Main sam- 
ples, and also the 2-degree Field Galaxy Redshift Survey 
(2dFGRS) data. In total, the sample comprises 893319 
galaxies over 9100 deg 2 . Baryon Acoustic Oscillations are 
observed in power spectra measured for different slices in 
redshift; this allows constraining the distance- redshift re- 
lation at multiple epochs. A distance measure at redshift 
z = 0.275 was achieved; but what is more important for 
our application of cosmography is the almost indepen- 
dent constraint on the ratio of distances, 



B = 



£>v(0.35) 
D v (0.2) 



= 1.736 ±0.065 



(33) 



which is consistent at l.lcr level with the best fit ACDM 
model obtained when combining the z = 0.275 derived 
distance constraint with the WMAP 5-year data. This 
measurement is particularly well suited for our purpose, 
because without it we could not apply BAO constraints 
on cosmographic approach. The usually derived distance 
measure from BAO are in the form r s (zd)/Dy, where 
fs (zd) is the comoving sound horizon at the baryon drag 
epoch, while Dy is defined as: 



Dv(zbao) = 



r BAO cdzV CZBAO 



H(z)J H(zbao) 



1/3 



(34) 



with H{z) the Hubble parameter. It is clearly obvious 
that the sound horizon quantity is incompatible with cos- 
mography, being it evaluated at at the drag epoch red- 
shift which is z w 1000, well out of the possible redshift 



applicability range of cosmography. On the contrary, the 
quantity Dy can be easily used for cosmography by sub- 
stituting for H(z) the appropriate series expression and 
evaluating it at z = 0.2 and 0.35, well inside the conver- 
gence redshift radius of cosmographic series. The quan- 
tity B can then be used for estimating cosmographic pa- 
rameters by minimizing the quantity 



Xbao(W) 



(6({0i)} - B obs f 



07- 



(35) 



V. COSMOGRAPHIC ANALYSIS 

As said in the previous sections, we have used an 
MCMC algorithm to perform the analysis of the multi- 
dimensional space of the cosmographic parameters. 

Our main purpose is to explore the possibility of ob- 
taining better constraints on the cosmographic parame- 
ters than offered by the literature. This is our motivation 
to use MCMC methods; as opposed to the approaches 
used in past examples in the same field, our choice offers 
the interesting possibility to introduce and manage easily 
possible physical constraints, thus making it possible to 
discover underlying relations and to realize a systematic 
report of the results. 



A. Preliminary discussion 

First of all, we try to answer a question: are there any 
physical constrains that can be imposed as priors to our 
fitting procedure? If yes: how do these constraints work 
on the estimation of parameters? 

The accuracy of the fitting parameters is of course a 
property derived from the likelihood function and it is 
quite independent of the method used for its maximiza- 
tion. The point that we want to investigate here is if there 
is any physical reasonable prior that can be introduce in 
such analysis for reducing the historically large uncer- 
tainties on cosmographic parameters (as described in the 
introduction). While a great part of these uncertainties 
surely depend on the series approach which cosmography 
is based on, we think that probably one could improve 
the analysis imposing some physical basic and general 
requirements. For example: does it make any sense to 
perform a minimization of x 2 without considering if the 
best fit parameters give well based physical information, 
i.e. without considering how physics is mapped in the 
parameter space (in our case: cosmographic parameter 
space)? 

This is a matter of great importance (and also their 
strong point in our idea) just while using statistical meth- 
ods like the MCMC: they are based on an algorithm that 
moves "randomly" in the parameter space: do all the 
points (even around the minimum x 2 location) tested 
by this algorithm satisfy some general physical require- 
ments, as a positive cIl or a positive H 2 1 These questions 
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would be of course useless if we were working with exact 
analytical expressions but, as log as cosmography is built 
on series expansions, we think it is necessary to analyze 
this aspect. 

As we will show below, even if we will be able to sen- 
sibly reduce uncertainties on cosmographic parameters, 
other problems arise (as, for example, the dependence on 
the maximum order series expansion, and on the data 
redshift range) that are intrinsic to the cosmographic ap- 
paratus. 

The main quantities we are going to fit are the Hubble 
parameter and the distance modulus; we are quite sure 
that these quantities will be physically well-based in the 
best fit location and around it. However, we have to take 
into account also quantities that arc not tested directly. 
For example, eventual constrains have to be given to the 
series expansions of H 2 (z) and d^z). The most general 
and obvious limitation is the positiveness of these two 
functions: for H 2 it is a natural mathematical property; 
and the same consideration is true for d^-, this being a 
distance. At the same time there is not a direct connec- 
tion between the positiveness of H and d^ ■ it is possible 
to have an H(z) function with varying signature, for ex- 
ample a sinusoidal one, and having a positive definite lu- 
minosity distance. We have also verified that, obviously, 
when fitting models in the redshift range of our data we 
automatically get a definite positive H(z), at least in the 
prescribed redshift range (thus excluding a bounce). 

Then another question arises naturally: does it make 
any sense to extend the priors (in our case the positive- 
ness of the physical functions) to any range of redshift 
well beyond the convergence radius of the series? The an- 
swer to this question is not so simple. Since we are work- 
ing with a truncated series expression we have to impose 
a "minimum requirement" for the positivity of our func- 
tions: it ranges only in the limited convergence redshift 
interval, < z < 1 or in the redshift range defined by the 
used data sets. For this reason we also perform a cosmo- 
graphic analysis both using the full (in a redshift sense) 
data samples and a cut sub-sample limited to z > 1. We 
also test what happens when extending positivity to all 
the possible redshift values, namely z > even if this 
would be rigorously right only if we were working with 
exact analytical functions, and to the maximum redshift 
value coming from any observational sample. 
Instead, in the y- redshift, things arc different: this quan- 
tity spans all the physical distances we are interested in, 
and it remains always inside the series convergence ra- 
dius, as < y < 1. So we can impose positivity on all 
the y-redshift interval without having forced final results 
and loss of information due to this choice. Anyway we 
also test results when the y-redshift is constrained to the 
range corresponding to the relative z-redshift range for 
comparing results. 

We have also found it is possible to add another con- 
straint: the requirement that VL m be positive and smaller 
than unity. Using the procedure we will describe in 
next section, where we will compare our cosmographic 



analysis with dark energy models described by the CPL 
parametrization, it turns out that m is a function of the 
cosmographic parameters. In the two dimensional case 
we have 



^m(<?o, jo) 



2 (jo - g0 - 2q 2 ) 
1 + 2j - 6g 



(36) 



while in the three dimensional case the expression be- 
comes quite longer, so we omit it for the sake of simplic- 
ity. Clearly, the above mentioned constraint on Q m can 
be enforced while running the MCMC algorithm. 
At the end, our priors are: 

• d L (z) > ; 

• H 2 (z) > ; 

• < fl m < 1 , 

applied on both z and y redshift ranges. 

Finally, for fully understanding and comparing our re- 
sults with literature, it is order now to underline the 
different use we have made of the jerk as compared to 
previous references. In [l[ and the quantity jo + Oo 
was estimated, where 



On = 1 



kc 2 

TT2 r ,2 

H a 



(37) 



As stated in § (|VII[) . in our case the spatial curvature 
is assumed to be k = so that we can derive jo values 
also from other different works. Of course when compar- 
ing our results to these ones, we have to consider that 
in those cases error bars correspond to the composite 
quantity jo + f2o- But we are confident that, consider- 
ing independent estimations of the curvature parameter 
and its related error, the largest contribution to errors is 
mainly attributable to the jerk parameter. On the other 
side, in Q, [38| and [39| the jerk parameter has the same 
definition as we use. Finally, in [1(| GRBs were the only 
data used (instead of SNela) . 



B. z-redshift 

Finally, if we consider the analysis with the common 
definition of z-redshift, looking at Tables (jH - (HJ) we can 
summarize the main results as follows: 

• with two cosmographic parameters, working with 
cut subsamples, namely considering only data with 
z < 1, always gives better results than using the 
complete sample and results are mostly uniforme. 
This is a direct consequence of the SNela sample 
having the largest weight in the chi-square min- 
imization because of the highest number of data 
points with respect to other data sets; 

• with three cosmographic parameters the differences 
between using total and cut samples decrease. Best 
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results always come from the cut samples but this 
time, differently from the previous point, they are 
related to the high-redshift extended priors, while 
assuming priors in the series convergence radius 
range does not satisfy the MCMC convergence re- 
quirements. This is somewhat expected from dis- 
cussion in § and possibly due to the fact that 
adding a third parameter in the restricted redshift 
range does not give any further improvement, while 
needing it for comparing cosmographic series with 
larger redshift ranges; 

• adding GRBs to SNela (it was interesting to ver- 
ify what happens when joining these data sets be- 
cause they constitute an homogeneous sample as 
they are both based on the definition of the distance 
modulus) when considering two cosmographic pa- 
rameters does not improve the analysis resulting 
in a very high chi-square value. When working 
with three cosmographic parameters the only sta- 
tistically good results come from the cut samples. 
But in these cases the weight of GRBs is low (only 
19 GRBs have redshift less than one) and no sta- 
tistically significative changes are detected in the 
estimation of cosmographic parameters; 

• the contribution of BAO data is almost statistically 
insignificant; 

• adding only Hubble data to SNela in the two pa- 
rameters case, while considering the total sample 
case, slightly increases chi-square value and pro- 
duces different values for cosmographic parame- 
ters. Introducing the prior on positiveness of phys- 
ical quantities on the Hubble data redshift range 
makes cosmographic estimations more similar to 
the SNcIa-only case. On the other side, when work- 
ing with three parameters we detect an higher influ- 
ence of Hubble data in the chi-square balance, with 
even largest (and previously undetected) deviations 
when applying the prior on positiveness z < 1.76. 
A further discussion on the reconstructed cosmo- 
logical scenario is given in § (|VT|) . 

With these considerations in mind, we can take a look 
to the cosmographic parameters coming from using only 
SNela and considering just two parameters: the decel- 
eration and the jerk. This is the most popular cosmo- 
graphic parametrization in the literature so that our re- 
sults and the advantages of the MCMC method can be 
easily judged and tested. We also underline that all the 
results in the literature are taken from sample of SNela 
where no redshift cut is applied. 

As we have pointed before, there is a dependence of re- 
sults from priors and from using cut or total sample which 
gets clear when looking to Fig. (JSJ) - If we move from 
the total sample and the prior applied on the total red- 
shift range z > (top panel), to the total sample with 
prior applied on the SNela redshift range (middle panel) , 
and finally to the cut sample with prior in the z-rcdshift 



convergence series radius, we detect a trend in the cos- 
mographic parameters: the deceleration starts far from 
usual values in literature (and ACDM one too) and moves 
towards them, while jerk starts from a value very near to 
the ACDM one and which agrees perfectly with most of 
the results in literature, and moves to smaller values. 

Considering the best chi-square values of the decelera- 
tion parameter, we remember that we are just using the 
Union2 data set from [30l |. which mixes many indepen- 
dent SNela observations so we may expect some differ- 
ences with results like the ones in [l[ . It is also interesting 
to underline that we agree with the analysis from [38| . 
which uses a completely physically different data set, the 
X-ray observations of galaxy clusters, together with Gold 
and Legacy supernovae. Their method is also different: 
while in [1[ they use the same series expansion approach 
as we do, in (381 ] an alternative numerical way leading 
to analytically exact functions was defined. The main 
difference in our work with respect to [l| is the use of 
physical constrains. Givin g ou r agreement with the exact 
analytical approach from [381 ] . we may argue that these 
constrains are an irrefutable requirement when working 
with cosmography as a series expansion tool. Finally, we 
perfectly agree at la confidence level with the results in 
Q, while only slightly agree with results in fioj coming 
from the application of cosmography to high redshifts 
using gamma ray burst data. 

About the jerk parameter it can be verified there is a 
good agreement with the analysis from X-ray emission 
from clusters of galaxies described in [38j and with the 
value coming from the Gold data set in |l|, j us t as in the 
deceleration case. 

At the same time we underline that it is also nude-eye 
visible that the error bars coming out from our analysis 
are actually narrower than in any other case (except for 
who derives cosmographic parameters from theoret- 
ical assumptions without using any observational data). 
If we consider that we are working with a physically ho- 
mogeneous sample (it is made of SNela only) but built 
from different (obscrvationally and technically reduced) 
data sets, this positive result may be strongly related to 
the chosen statistical method, thus MCMC are showing 
their good nature with respect to other popular fitting 
procedures. 

When moving to the three dimensional analysis, in- 
cluding the possibility of constraining the snap param- 
eter, we find very different results. If we consider only 
the SNela sample, now the deceleration parameter is less 
negative and the jerk even changes sign becoming nega- 
tive. We find a very weak match with past values from 
literature mainly due to their large errors. On the con- 
trary, using MCMC makes us able to narrow down the 
confidence levels strongly. 

The discussion cannot be complete without considering 
the value of the snap. Unfortunately we have not many 
examples in the literature we can compare with: from |39j 
and [1 01 ] we have some strongly positive values with per- 
centage errors much larger than the 200%, while our anal- 
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ysis gives a negative value with narrower constraints. We 
remind the reader that the snap value in ACDM should 
be equal to —0.35; in this sense, the sign we have found 
shows a clear trend in this parameter and in the cosmo- 
logical properties subtended by it. 

When adding Hubble data, we have two different sit- 
uations related to the different applied priors: when the 
positiveness is for z > we find the main differences for 
the jerk, which is much more negative, and for the snap, 
less negative; on the other side, when the positiveness if 
for z < 1.76, differences arc more pronounced: the decel- 
eration is very near to —1, while the jerk and the snap 
are positive and larger than before, 3.134 and 4.399 re- 
spectively. It is not an easy task to understand why this 
happens; moreover no correspondence with past values 
in literature is found. 



C. y-redshift 

If we work with the y-redshift defined in [l[ to solve 
the convergence problem of the series, looking at Ta- 
bles (|III[) - (|IV[) we can find some interesting changes: 

• as expected from the definition of the y-redshift, 
differences both among total and cut samples and 
using different y-rcdshift prior ranges are strongly 
softened because we are always well inside the con- 
vergence radius of the series; 

• results are always consistent and uniform, showing 
only small statistical variations, when considering 
two and three cosmographic parameters separately, 
but cosmographic estimations are different between 
the two and the three dimensional approach; 

• adding high order y-rcdshift data (GRBs) to SNcIa 
does not change results in a significant statistical 
way. We have to remind that the y-redshift scries 
show problem in reconstructing the distance modu- 
lus for y w 0.4 -j- 0.5, which is well under the SNcIa 
maximum rcdshift even. So, we feel that this fea- 
ture mostly depends on the intrinsic impossibility 
of y-rcdshift scries to match the analytically exact 
physical behavior of distance modulus; 

• BAO has no weight on the estimation of cosmo- 
graphic parameters; 

• as happened in the z-rcdshift case, here too Hub- 
ble data are the only exception in the uniformity 
of results exhibiting a stronger weight in the cos- 
mographic estimations with respect of GRBs and 
BAO. While the contribution of the cut Hubble 
sample is not significant, the total one, extended 
up to z = 1.76, i.e. y = 0.638, produces the 
largest deviation. As the Hubble sample does not 
extend the y-rcdshift range in a much sensible way 
for explaining such a deviation, we can only argue 
that some peculiarities in Hubble high redshift data 



are present also inducing the difficulties in fitting 
them as discussed in § (|VI[) . A corroborating test 
for this explanation is the not reached convergence 
for MCMC in the three dimensional analysis when 
the total Hubble data sample is considered. 

When comparing our results with the literature, the 
only previous examples where the y-redshift was used 
can be found in Q. As in our case, when moving from 
z-redshift to y-redshift they found an increase in errors 
(which are bigger than ours, anyhow) of the cosmo- 
graphic parameters, more consistent in the case of the 
jerk, as we have verified. This could be given by the sub- 
stantial down-stretching of the redshift data when us- 
ing this different parametrization; even if the physical 
content of the data is not altered, it is more difficult 
to extract cosmographic results when the data are flat- 
tened on the new redshift axes. The z-redshift range of 
the Union2 data set, < z < 1.4, goes over to the y- 
rcdshift interval < y < 0.584; Hubble data convert 
from < z < 1.76 to < y < 0.638, and GRBs from 
< z < 6.6 to < y < 0.871. It is clear that if we con- 
sider the great contribution from SNela to the chi-square 
function, we are fitting our physical distances on a limited 
range (low y-redshift) with respect to the total conver- 
gence radius (y = 1) so that higher order parameters (as 
the jerk and the snap) can be not as well estimated as 
lower order ones. Thus while for the deceleration param- 
eter errors still remain quite centrally based, for the jerk 
case we have a long tail for high positive values (nearer 
the ACDM value but only matching it at 2a confidence 
level) while for the snap we have a tail for large negative 
values matching the corresponding ACDM value. This 
may be considered a trend which one eventually needs to 
test and confirm if having higher order data. 

The agreement is good for most of the past works we 
have considered, with the greatest number of successful 
matchings for the deceleration parameter than for the 
jerk. Anyway, the best results at all, come from the case 
where we use the full data samples and applied the prior 
in the range depicted by them. In this case also the decel- 
eration and the jerk match well with the corresponding 
ACDM values. 



VI. EQUIVALENT CPL MODEL 

To round up results from simulations on cosmographic 
parameters, we can go deeper into the analysis and try 
to our results with a specific cosmological model. It is 
easy to check the link between them. In this paper we 
have used the CPL parametrization for dark energy given 
by Eq. (|17|) which depends only on three parameters, 
(Q m , w Q and w a ). Using Eqs. - - (ED we can 
relate cosmographic parameters and CPL ones. We have 
different possibilities depending on the number of cos- 
mographic parameters we are working with and on how 
many CPL parameters we are going to consider as free. 
In the end we have: 
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• with two cosmographic parameters, (go, jo), we can 
derive news about a Constant Dark Energy (CDE) 
model (i.e. w a = 0), with: 




Figure 4. Hubble parameter H(z) and cosmographic pa- 
rameters q(z), j(z), s(z) in different cosmological scenarios: 
ACDM model (solid gray line); Komatsu et al. (2010) model 
(dashed gray line); Perivolaropoulos et al. (2008) model (dot- 
ted gray line); Model I (long-dashed black line); Model II 
(short-dashed black line); Model III (dotdashed black line); 
Model IV (dotted black line). 



to m (q ,j ) 

w (q ,jo) 



2(j -g -2g 2 ) 

1 + 2j - 690 
1 + 2 j - 6g 
-3 + 6g 



(38) 



• with two cosmographic parameters, (go, jo), we can 
derive news also about a Dynamical Dark Energy 
(DDE) model, (i.e. w a ^ 0), leaving il m free, with: 



wo (go, jo, ^m) 



l-2g 



3(-i + n m ) ' 

Wa(qo, 30, torn) = ot TTTT^ ( _2 io(-l + to m ) - 2g 

3(-i + \i m y 

x (l + 2g ) + a„(-l + 6g )) ; (39) 

• with three cosmographic parameters, (go,jo,So), 
we can derive news about a Dynamical Dark En- 
ergy (DDE) model, with fl m depending on cos- 
mographic parameters, i.e. £l m = Cl m (qo, jo, So)- 
The same holds true for wq = wo(go, jo, so) an< ^ 
w a = w a (qo 1 jo, so). These expressions are much 
longer than previous ones and we do not write them 
here for the sake of brevity. 

The most important thing is that with these relations all 
the statistical properties of these parameters (median, 
error bars, etc.) can be easily extracted from the cos- 
mographic samples we have obtained from MCMC sim- 
ulations. Final results of this dark energy analysis are 
reported in the tables on previous pages. 

The first main result concerns Hubble data and their 
relation with the CDE model available from using cos- 
mographic series with two parameters: when using z- 
redshift, this dynamical model is always unable to fit 
Hubble data. Even if we consider the model with param- 
eters directly obtained from using this kind of data, we 
are only able to fit few points at very low redshift and 
not all the sample. Moreover, when moving to the CDE 
models with y-redshift and to the DDE models derived 
from three cosmographic parameters using Hubble data, 
we can easily verify that the derived cosmological param- 
eters (and the cosmographic one too, as we pointed out 
in §. (jVj) are always very different from the ones derived 
from other kinds of observational data. This feature can 
be interpreted in two different ways: 

• it depends on the procedure which Hubble data are 
obtained from; 

• there is an intrinsic information about the dynam- 
ical model in Hubble data. 

Surely, for the two cosmographic parameters bet- 
ter way of proceeding may be to obtain Q m from alter- 
native ways and then join results with those one from 
cosmography as we discuss later in this section. 
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If wc consider the CDE model derived from two cos- 
mographic parameters and the DDE models derived from 
three cosmographic parameters with y-redshift, we can 
see that: 

• fl m is lower than usual, being « 0.07, and matching 
at the top lcr confidence level the usual estimation 
of f2 m ~ 0.16 obtained considering plasma contri- 
bution; 

• if we also add Hubble data, then the value for the 
matter content grows up to « 0.40 in the CDE 
case and to ss 0.36 in the DDE case. In both cases 
the bottom lcr confidence limit matches the mat- 
ter content value (plasma included) and the total 
ACDM one. 

It is much intriguing to find such small values. Gen- 
erally from SNela one is able to extract a slig htly over- 
estimated value of fl m ~ 0.30 in [4(j; in unifying 
SNela, BAO and WMAP7, a value for fi m ~ 0.245 was 



obtained; while in [38( it is fi m = 0.306_ 'o40- We under- 
line that any conclusion about the real possibility of such 
universe dynamics is out of the purpose of this paper. 
We have to consider that we are always working with an 
approximated approach (series expansion) and that this 
conclusion should be corroborated by a more extended 
analysis of all the other dynamical aspects of universe 
(formation of any gravitational structure, explanation of 
acceleration of universe) . Then we need to refine the use 
and the approach of cosmography (what we are intended 
to do with a series of papers). 

In the same cases as before, for what concerns the dark 
energy parameters, we have: 

• for the CDE models, Wo ~ —0.82, excluding the 
phantom line crossing at la level; while adding 
Hubble data it becomes much more negative, wq ~ 
-1.59; 

• for the DDE models, wo « —0.715 and w a « 0.01 
(~ —0.16 at la error); adding Hubble data we get 
w a « -1.272 and w a w -0.326. 

The DDE models from two cosmographic parameters 
have £l m fixed at the value obtained in (ill ], while the 
dark energy parameters are: 

• with the z-redshift, wq w —1 and w a s» —0.71; 
adding Hubble data they change to ido ~ —1-4 and 
w a ~ —1-2, but they also have high chi-square val- 
ues. About the w a parameter, the question is very 
subtle and tricky: we have a negative value for it 
very far from the usual theoretical expected values 
and almost of the same order of wq . 

• with the y-redshift, wo ~ —1 and w a « —0.56 
which match well the results from [ll|, wo = —0.93 
and w a = —0.41. Adding Hubble data they change 
to wo ~ —1.25 and w a ~ 1.2. In this case, on the 
contrary, the value of w a is of the same order of wq 



but positive so that wq + w a ~ and agrees with 
results from [4l[ where starting from Union SNela 
sample find wq = —1.4 and w a = 2. 

To test all this results we can also verify that consider- 
ing the evolution of Hubble parameter or of the deceler- 
ation parameter, we do not find particular problems, but 
rather some intriguing properties. Using Eqs. (fT8|) we can 
derive the functions q(z,Cl m ,Wo,w a ), j(z,Q m ,Wo,w a ) 
and s(z, Q TO , wq, w a ). In this way, giving the cosmolog- 
ical model obtained by our cosmographic analysis and 
parameterized by the CPL prescription, we can infer the 
redshift dependence of the cosmographic parameters. In 
Fig. ((3]) we have plotted the analytical behavior of H(z), 
q(z), j(z) and s(z) for five different cosmological models: 



ACDM: ft m = 0.30, 



w 



-1, 



0; 



• best cosmological model from [15J: Q r , 
w = -0.93, w a = -0.41; 



• best cosmological model from [41|: f2 r 

w = -1.4, w a = 2; 



0.245, 



0.30, 



• Model I, CDE model from two dimensional cosmo- 
graphic analysis, in y-redshift applying priors only 
in the redshift range depicted by using all the data 
samples: Q m = 0.399, wo = —1.581; 

• Model II, DDE model from two dimensional cosmo- 
graphic analysis, in y-redshift applying priors only 
in the redshift range depicted by using all the data 
samples: Q m = 0.245, w Q = -1.256, w a = 1.220; 

• Model III, DDE model from three dimensional cos- 
mographic analysis, in z-redshift applying priors in 
the full z-rcdshift range and using all the data sam- 
ples: fi m = 0.151, w = -0.652, w a = -0.207; 

• Model IV, DDE model from three dimensional cos- 
mographic analysis, in y-rcdshift applying priors 
only in the redshift range depicted by using all 
the data samples: fi m = 0.364, w = —1.272, 
w a = -0.326; 

The first thing we can say is that our models can be 
different from the ACDM scenario for low redshifts, but 
they are equivalent to it for high redshifts. This is not 
the case of the model in [4l|, which is completely dif- 
ferent from ACDM model at any redshift. Of course, in 
our case, we have the pleasant property of that some 
properties at high redshift (such as CMB theory and 
data) remain untouched; but, of course, the necessity 
of re-discussing what happens at low redshifts cannot be 
waived. 

Considering Model I, for example, it has a transition 
from positive to negative q(z) at z = 0.44 while this 
happens at z = 0.68 in a ACDM scenario. At the same 
time, transition from a matter dominated epoch to a dark 
energy dominated one happens at z ~ 0.1 in our model, 
while it is at z = 0.32 in ACDM. 
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Considering Model II, for example, it has a transition 
from positive to negative q(z) at z = 0.62 and from dark 
matter to dark energy domination at z ~ 0.42. Even 
if it has dynamical parameters different from the ACDM 
ones, it also exhibits gross properties that are very similar 
to the consensus model, but at the same time the jerk and 
the snap behavior are much different only converging to 
ACDM values at high redshifts. 

Model III and IV arc different from previous models 
because dark energy fluid (given our CPL parameters val- 
ues) has always been an important ingredient of universe 
dynamics with negative pressure (w a is negative). The 
Model III seems to be the most problematic showing a 
transition from positive to negative q(z) at z = 1.52 and 
from dark matter to dark energy domination at z ~ 1.23. 
On the contrary, Model IV shows more normal prop- 
erties, as a transition from positive to negative q(z) at 
z = 0.53 and from dark matter to dark energy domina- 
tion at z ~ 0.16 also exhibiting very similar trends for 
the jerk and the snap. 



VII. CONCLUSIONS 

In this paper we are intended to perform a very de- 
tailed analysis of Cosmography, to state its degree of us- 
ability and cxtendibility when used as a model to universe 
dynamics. 

First of all we have shown that results from a cosmo- 
graphic approach are sub-judice because of the observa- 
tional errors of cosmological data one works with. When 
working with cosmography one has to choose the highest 
order of expansion for the series of the a(t) function so 
stating how many cosmographic (i.e. cosmological) pa- 
rameters it is possible to extract from the analysis. But 
we have seen that differences between using a two, three 
or four dimensional scries are well inside the errors given 
by data, so that (at the moment) any possibility of ex- 
actly stating what kind of scries is better to be used is 
out of question. 

From statistical considerations, we have seen that 
working with the series expansion approach, it is pos- 
sible to derive up to a maximum of three cosmographic 
parameters, the deceleration, the jerk and the snap. Fi- 
nally, the statistical method we have chosen to work with, 



MCMC, gives us the possibility to sensibly narrow down 
the confidence errors on cosmographic parameters with 
respect to other results in the literature (coming from the 
same series approach). 

Then MCMCs give also the possibility to manage pri- 
ors on physical quantities. Comparing results with the 
literature we see that inside the narrower errors our esti- 
mations still agree with many works; what is more inter- 
esting, is that they also agree with a different approach 
leading to numerical evaluation of an exact analytical 
function of distance modulus ([38j)- 

What is even more interesting is what seems to come 
out from the analysis of dark energy parameters (in the 
CPL parametrization) related to the set of cosmographic 
parameters we have found out. We also know that Cos- 
mography could be used as an independent tool to detect 
possible deviations from the ACDM model; but until now 
nobody has checked what the obtained cosmographic sets 
means in terms of dynamical information. 

We have tried to give an answer, but this is still quite 
puzzling: cosmographic sets can be related to DE mod- 
els which are quite slightly different from ACDM (at least 
at low redshifts). At the same time, smaller differences 
in the Cosmographic parameters values can give great 
differences in the dynamical cosmological model. For ex- 
ample we have seen that in the case of two cosmographic 
parameters, when working with only SN, we have a good 
agreement with literature. But these values correspond 
to a constant dark energy model with a low matter con- 
tent. At the same time, this model is not able to fit 
Hubble data at high rcdshift. 

Finally, very different global dynamics properties, can 
produce the same gross trends for Hubble and deceler- 
ation parameters (as much the jerk and the snap ones) 
pose the question of a possible cosmological dynamical 
degeneracy in the universe, which also makes more diffi- 
cult to eventually detect alternative scenarios to ACDM. 

In this optics developing cosmography is most impor- 
tant giving the possibility to discriminate among models 
using absolute values of cosmographic parameters. We 
think that all these questions could be solved improving 
the series expansion approach of Cosmography: we are 
studying the possibility to define an Exact Cosmography 
which could be solve the cosmographic degeneracy. A 
detailed description of this approach will be given in a 
forthcoming paper. 
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Appendix A: Cosmological distances in 
Cosmography 

The physical distance traveled by a photon emitted at 
time and absorbed at the current epoch to is 



D 



dt = c(t — t*), 



(Al) 



so that inserting t* = to — D/c into Eq. gives us an 
expression for the rcdshift as a function of to and D/c, 
i.e. z = z(D). Solving with respect to D, up to the fifth 



order in z, gives us the desired expansion for D(z): 
D(z) = ^{V° + Vlz + Vlzi + Vlz 3 
+Vtz A 



0[z 5 ]} 

In the latter we have defined 
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The Taylor series expansion of the quantity D(z) is a 
building block for other quantities to be derived in the 
following sections as required by the cosmographic ap- 
proach. 

In typical applications, one is not interested in the 
physical distance D(z), but in other definitions. In our 
case, the luminosity distance, 



D, 



a (to 



-(o(to)r (I>)), 



a(to - D/c) 

has a particular relevance. In Eq. (|A3[) we have used 
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Using Eq. (|A2[) . and after some lengthy algebra, we ob- 
tain 
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One now must rewrite tq(D) as a function of z using 
Eq. (|A2[) . Inserting the result into Eq. (|A3|) . one obtains 
the fifth order approximation for the Hubble free lumi- 
nosity distance dh = Dl(z)/(c/Hq) as a function of the 
redshift z: 



d L {z) = Vlz + Vl 



V 2 z 3 



vi, + 



V\ z 5 . 
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In the latter we are using the definitions 
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The distance modulus is 

tiv) = r^r ( lo § y + Ml y + M "y 2 + M3 y 3 + M V) , 
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The Hubble parameter is 

H(y) =Ho-(Ul + Uly + H 2 y 2 + H 3 V 3 

+ n 4 Q y 4 ) (B5) 

where 



Ho — 1, 



H Q = l + q , 



(B6) 
(B7) 



Previous relations in this section have been derived 
for any value of the curvature parameter; in (42|, it is 
shown that ranging the curvature parameter in the in- 
terval [—1,1] has negligible effects on the estimation of 
cosmographic parameters, so that, in the following of the 
paper, we will assume a flat universe and use the simpli- 
fied versions with k = 0. 



Appendix B: Y redshift relations 

We write here all the needed relations expressed in 
term of the y-redshift for k = 0. The Hubble free lumi- 
nosity distance <Zl(j/) is: 

d L (y) = V\y + Vly 2 + V\ y 3 + V\ y 4 + V\ y 5 , (Bl) 

where 

vi = i, 

»i = |[i-«b], 
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The square of the Hubble parameter is 

H 2 (y) = H 2 (Hi' + Hl x y + H 2 Q 2 y 2 + H 2 ' 3 y 3 

+ n 2 Q A y 4 ) , (Bii) 
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Appendix C: H derivatives 
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Here we report convenient expressions for the conver- 
sion from higher order derivatives in time to higher order 
derivatives in redshift: 



Appendix D: Square Hubble parameter 
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The expansion series of the square Hubble parameter, 
H 2 , will be a useful tool as well. It is a simple and just 
a long matter of algebra to compute its derivatives with 
the same procedure as above. The final results are 
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Table I. Results in the two dimensional parameter space (50, jo) and z-redshift. Column (1): simulation identification in the 
text with the used data (S = SNela, G = GRBs, H = Hubble, B — BAO; the index (c) refers to the redshift-cut sub-sample, 
with z < 1) and the applied priors (we indicate the redshift range where positiveness of quantities described in § ()V A|) is 
applied). Columns (2) - (3) - (4): best fit estimations of cosmographic parameters with la confidence level. Column (5): chi 
square and reduced chi square values. Columns (6) - (7) - (8) - (9): best fit estimations of dark energy (CPL) parameters with 
la confidence level for the constant dark energy model and the dynamical one with fixed fi m as described in § (I VII) . 



Data 



Best fit parameters 



CDE 



DDE (fi m = 0.245) 



Jo 



X 2 (xVd.o.f. 



n, 



Wq 



w 



z>0 
z < 1.4 

z > (c) 
z < 1^ 



-1.014 



+0.008 
-0.013 



-0.738±g:gg 

-i.oi7ig;S? 

-0.669l° : °g 



+0.063 
-0.035 

-0.099 
-0.085 

-0.082 
-0.045 



0.284t°;^ 



1.060 
0.372; 
1.072 



622.11 (1.123 
573.95 (1.036 
599.21 (1.122 
529.97 (0.992 



0.004 



+0.005 
-0.003 



■1.013 



+0.008 
-0.014 



-1.337±g;S?I -1-295 



no5 +0 009 

U.UUO_ 004 

n05 +0 - 006 

U.UUO_ 003 

01 9+ - 027 
u.u±y_ 014 



-0.831_°;°__ 

-i.oieiroi? 

-0.796lg;gS 



, nq o + 0.039 
i.u»o_ 038 

1 qqo+0.009 

-l.OO3_g 014 

1 n^+ - 043 
-i.uoo_o 046 



+ 0.005 
-0.011 



D eeo + 0.059 

1 OQc;+0.005 
-i.z,ao_ on 

D 790 + 0. 063 
-U. I zo_ 064 



S-G 



Jo 



X 2 (x 2 /d.o.f. 



n, 



Wq 



Wo 



z>0 
z < 6.7 

z > (c) 
z < 1^ 



-1 004+ 003 

l.UU1_o.oo5 

-0 ^n7 +0 026 
v.oui _ .o3i 

-1 016+ 009 

l.UlD_o.oi7 

-0 6 73+ 045 

U.O(0_o.o50 



1.016 



-0.016 
-0.011 



010+ ' 034 

U.U±U_ .o25 



1.068 



-0.084 
-0.042 



292+ - 149 



1166.07(1.872 
1015.14 (1.629 
602.03 (1.089 
533.20 (0.964 



0004+ 0038 

U.UUU^_ .0032 

001+ 002 
u.uui_ 001 

04+ - 007 

U.UU<±_o.oo3 

20+ - 027 

u - uzu -0.015 



-1 no3 +0 002 

l.UUO_o 003 

-0 672+ 017 

U.U/ ^_o.020 

1 m k+o.009 

-l.U10_o.oig 

-0 7QQ+0-035 
U. / »»_ .049 



1 090 + O.OO3 
-l.OZO_ 004 

-0.889_K 

1 000+0.008 
-1.ooo_0.01g 

-1 nS6 +0 - 040 

l.UOU_ 044 



1 9QO + 0.014 
-l.^»O_ 023 

fl K7C + 0.030 
-U.OI O_o.o36 

1 9Qc;+0.005 
-i.^»o_ 010 

O 797+0.065 
"U. z _q.062 



S-H 



qo 



Jo 



X 2 (x 2 /d.o.f. 



Wq 



WO 



z>0 
z < 1.76 

z > (c) 
z < 1^ 



-1 15Q+ 0059 

l-lO»_0.048 

-0 fiSfi+ 052 

U.UOU_o 056 



■1.016 



+0.010 

-0.016 



-0 667+ 048 
u.uu* _ .o5i 



1.859_o; 2 3y 

91 q+ - 152 

U.^l»_o 095 

1 069+ - 081 
i.uoa_ 04i 

u.288_o.io3 



717.88 (1.268 
637.70(1.127 
606.06(1.118 
532.80 (0.983 



o.o56±g : gi| 

o.oi7ir ?i 

U.UU1_o.oo3 

+0.028 
-0.016 



0.022 



1 170+ 063 
-1.1 /U_ .o53 

ft 771 +0.040 
-U. I I l_o.052 

l.Uld_o.oi 7 

-0 7Q7+0.039 
u - ' vl -0.049 



-1 4R4+ - 052 

l.Wi_ 042 

-1 no3 +0 - 046 

l.UUO_ 049 

1 oqn+0.008 
-i.oo»_g 014 

, qoi+0.042 
l.U01_Q 046 



1 9CO + 0.010 
-l.ZOO_ on 

-0.689t° ;°l 

1 9Qc;+0.005 
-i.z»o_0 010 

-0 714+ 065 

U. / 11_o.062 



S-G-H 



<lo 



Jo 



X 2 (x 2 /d.o.f. 



n, 



Wq 



Wo 



z>0 
z < 6.7 

z > (c) 
z < 1^ 



-1.004 
-0.476 



+0.002 

-0.004 

+0.010 
-0.018 



-1 017+ - 010 

1.U1 1 _ .017 

-0 6 67+ 050 

U.ULH _ .052 



1 mc + 0.015 
l.Uld_ .009 

-0 019+ 017 
u.uiy_ .oii 

1.074tr o f 6 
0.285t 9 ;^ 2 



1267.78 (1.996 
1050.51 (1.654 
609.44 (1.086 
535.94 (0.955 



0.0005 



+0.0008 
-0.0004 



-1.003 



+0.002 

-0.003 



-1.328 



+0.002 

-0.004 



-1.293 



+ 0.003 
-0.006 



001 +0 002 
u.uui_o.ooi 

005+ - 007 

U.UUO_o.o04 



0.021 



+0.030 
-0.015 



-0 652+ 007 
u.uoz_ 012 

-1 016+ 0010 

l.UlU_ .oi8 

-0 7 96+ 040 

U./»U_ .o52 



-o.862±S:S_S 

-1.340t° ;° 1 9 
1 rm+ - 044 

l.UOl_ 046 



-0 543+ - 011 
u.o^o_ 020 

1 9Qc;+0.005 
-i.zao_0.0io 

-0 714+0-068 
U. / 1<±_o.o64 



S-G-H-B 



qo 



Jo 



X 2 (x 2 /d.o.f. 



Wq 



Wq 



z>0 
z < 6.7 

z > (c) 
z < 1^ 



-1 004+ 002 

l.UU^_o.oo3 

-0 477+0.010 

u -^' ' -0.018 



1 01 4+ 0016 

l.U14_ 008 



+0.017 
-0.011 



■1.016 



+0.009 
-0.016 



-0.018 
1.072 



-0 6 67+ 049 

U.VVI _o.Q50 



-0.080 
-0.043 



0.282; 9 ;^ 3 



1267.61 (1.993 
1051.21 (1.653 
609.08 (1.084 
536.17(0.954 



0005+ 0008 

U.UUUO_ .0004 

001+ - 002 
u.uui_o.ooi 

005+ - 007 

U.UUO_ .003 

20+ - 028 
u.uzu_ .oi5 



-1 n03 +0 002 

i.uuo_ 003 
u.uoo_ 012 
-1 016+ 009 

l.UlD_ .oi7 

-0 7Q5+ 040 



1 090+0.002 
-i.ozo_ 003 

-0.863tr o ? 9 6 

1 qqn+0.008 
-i.oo»_Q 014 

-1 n30+ - 044 

l.UOU_ 044 



1 9QO + 0.003 
-l.z»o_ 006 

-0 544+ 011 

U.O^_o.o20 

-1 9q4+ 004 
- L -^ a ^-o.oio 

n 71 c+0.065 
-U. ( lO_ .063 



S-H-B 



qo 



Jo 



X 2 (x 2 /d.o.f. 



Wq 



Wq 



z>0 
z < 1.76 

z > (c) 
z < 1^ 



1 1 C9+0-052 
-l.lJZ_o 050 

-0 6 37+ 050 

U.UOf _ .059 



1 oiQ+0.308 
1-Ol»_o 294 

D 990+0. 163 
U.zzo_ g2 



1.017 



+0.010 
-0.017 



1.074- 



-0.083 
-0.045 



-0 6 69 +0 048 
u.uuy_ .05i 



0.288t 9 ;^ 2 



718.55(1.267 
638.01 (1.125 
606.23(1.116 
532.98 (0.981 



0.055 



+0.015 
-0.017 



o.oisto 1 ;^ 



+0.007 
-0.004 



0.005 
0.021 



+0.028 
-0.015 



-1 165+ 057 

1.1DO_ .055 

77Q+0.038 
"U. / / O_ 054 

-1 016+ 0010 

l.U10_ .oi7 

-0 707+0.040 

u - ' vl -0.048 



-1 459+ - 046 

l.U»_0 045 

-1 004+ ' 044 

l.UU £ ±_Q 051 
1 qqn+0.009 

1.003_g 014 

1 n^9+ - 042 

-l.UOZ_o 045 



1 9CO + 0.009 

-i.zoo_ 010 

-0 691+ 054 
u.oai_ .056 

1 9Qc;+0.005 
-i.z»o_ 010 

D 71O + 0.065 
-U.a8_o.o63 



19 



Table II. Results in the three dimensional parameter space (qo, jo, so) and z-redshift. Void rows correspond to MCMC which 
do not pass the convergence test. Column (1): simulation identification in the text with the used data (S = SNela, G = 
GRBs, H — Hubble, B = BAO; the index (c) refers to the redshift-cut sub-sample, with z < 1) and the applied priors (we 
indicate the redshift range where positiveness of quantities described in § (|V A|) is applied). Columns (2) - (3) - (4): best 
fit estimations of cosmographic parameters with la confidence level. Column (5): chi square and reduced chi square values. 
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Table III. Results in the two dimensional parameter space (qo,jo) and y-redshift. Column (1): simulation identification in the 
text with the used data (S = SNela, G = GRBs, H = Hubble, B — BAO; the index (c) refers to the redshift-cut sub-sample, 
with z < 1) and the applied priors (we indicate the redshift range where positiveness of quantities described in § ()V A|) is 
applied). Columns (2) - (3) - (4): best fit estimations of cosmographic parameters with la confidence level. Column (5): chi 
square and reduced chi square values. Columns (6) - (7) - (8) - (9): best fit estimations of dark energy (CPL) parameters with 
la confidence level for the constant dark energy model and the dynamical one with fixed fi m as described in § (IVI|) . 
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Table IV. Results in the three dimensional parameter space (qo, jo, so) and y-redshift. Void rows correspond to MCMC which 
do not pass the convergence test. Column (1): simulation identification in the text with the used data (S = SNela, G = 
GRBs, H — Hubble, B = BAO; the index (c) refers to the redshift-cut sub-sample, with z < 1) and the applied priors (we 
indicate the redshift range where positiveness of quantities described in § (|V A|) is applied). Columns (2) - (3) - (4): best 
fit estimations of cosmographic parameters with la confidence level. Column (5): chi square and reduced chi square values. 
Columns (6) - (7) - (8)): best fit estimations of dark energy (CPL) parameters with la confidence level for the dynamical dark 
energy model as described in § (|VI[) . 
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Figure 5. Comparing deceleration and jerk values from the literature with results from: (First panel.) SNela only, two 
dimensional cosmography and prior on z-redshift range z > 0; (Second panel.) SNela only, two dimensional cosmography and 
prior on z-redshift range z < 1.4; ( Third panel.) cut SNela only, two dimensional cosmography and prior on z-redshift range z < 
1; (Fourth panel.) cut SNela only, three dimensional cosmography and prior on z-redshift range z > 0. References in numerical 
order as they appear in the figure: (1) - John (2004); (2) - Rapetti (2006), Gold SN subsample; (3) - Rapetti (2006), Legacy 
SN subsample; (4) - Rapetti (2006), X-ray clusters subsample; (5) - Rapetti (2006), all subsamples; (6) - Poplawski (2006); 
(7) - Cattoen (2007), Gold SN in z-redshift; (8) - Cattoen (2007), Gold SN in y-redshift; (9) - Cattoen (2007), Legacy SN in 
z-redshift; (10) - Cattoen (2007), Legacy SN in y-redshift; (11) - Capozziello, Izzo (2008). Solid horizontal line is our best fit 
result; dotdashed horizontal lines show its la confidence level; dotted horizontal line shows ACDM expectation. 
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Figure 6. Comparing deceleration and jerk values from the literature with results from: (First panel.) all cut data, two 
dimensional cosmography and prior on z-redshift range z < 1; (Second panel.) SNela and Hubble data, three dimensional 
cosmography and prior on z-redshift range z > 0; (Third panel.) SNela and Hubble data, three dimensional cosmography 
and prior on z-redshift range z < 1.76. References in numerical order as the appear in the figure: (1) - John (2004); (2) - 
Rapetti (2006), Gold SN subsample; (3) - Rapetti (2006), Legacy SN subsample; (4) - Rapetti (2006), X-ray clusters subsample; 
(5) - Rapetti (2006), all subsamples; (6) - Poplawski (2006); (7) - Cattoen (2007), Gold SN in z-redshift; (8) - Cattoen (2007), 
Gold SN in y-redshift; (9) - Cattoen (2007), Legacy SN in z-redshift; (10) - Cattoen (2007), Legacy SN in y-redshift; (11) - 
Capozziello, Izzo (2008). Solid horizontal line is our best fit result; dotdashed horizontal lines show its la confidence level; 
dotted horizontal line shows ACDM expectation. 
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Figure 7. Comparing deceleration and jerk values from the literature with results from: (First panel.) SNela data, two 
dimensional cosmography and prior on y-redshift range < y < 1; (Second panel.) all data, three dimensional cosmography 
and prior on y-redshift range < y < 1; (Third panel.) all data, three dimensional cosmography and prior on y-redshift range 
< y < 0.871. References in numerical order as the appear in the figure: (1) - John (2004); (2) - Rapetti (2006), Gold SN 
subsample; (3) - Rapetti (2006), Legacy SN subsample; (4) - Rapetti (2006), X-ray clusters subsample; (5) - Rapetti (2006), all 
subsamples; (6) - Poplawski (2006); (7) - Cattoen (2007), Gold SN in z-redshift; (8) - Cattoen (2007), Gold SN in y-redshift; 
(9) - Cattoen (2007), Legacy SN in z-redshift; (10) - Cattoen (2007), Legacy SN in y-redshift; (11) - Capozziello, Izzo (2008). 
Solid horizontal line is our best fit result; dotdashed horizontal lines show its la confidence level; dotted horizontal line shows 
ACDM expectation. 



